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Abstract 

We investigate an infinitesimally thin cylindrical shell composed of counter-rotating dust parti- 
cles. This system was studied by Apostolatos and Thorne in terms of the C-energy for a bounded 



> 

o 

^sl ■ domain. In this paper, we reanalyze this system by evaluating the C-energy on the future null 



infinity. We find that some class of momentarily static and radiation-free initial data does not set- 
tle down into static, equilibrium configurations, and otherwise infinite amount of the gravitational 



^jj- radiation is emitted to the future null infinity. Our result implies the existence of an instability in 

this system. In the framework of the Newtonian gravity, a cylindrical shell composed of counter- 

cn ' rotating dust particles can be in a steady state with oscillation by the gravitational attraction 

^' 

CN ' and centrifugal repulsion, and hence a static state is not necessarily realized as a final state. By 



contrast, in the framework of general relativity, the steady oscillating state will be impossible since 



f~~- I the gravitational radiation will carry the energy of the oscillation to infinity. Thus, this instability 

O 

has no counterpart in the Newtonian gravity. 
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I. INTRODUCTION 

In order to understand the relativistic gravitational phenomena, the existence of isome- 
tries is often or usually assumed, since the Einstein equations are very complicated system of 
quasi-linear partial differential equations. The simplest but useful assumption in relativistic 
astrophysical situations is that of spherical symmetry. A shortcoming of this assumption is 
that there is no freedom of gravitational radiation. It is hardly possible to extract any effects 
of gravitational radiation within spherically symmetric systems. In contrast, although the 
cylindrical symmetry might rarely appear in relativistic astrophysical situations, it has a 
degree of freedom of gravitational radiation known as Einstein-Rosen gravitational waves 
PJ. This system has been studied in connection with the gravitational waves [2, y, Ij]- 

The cylindrically symmetric system has been investigated also in connection with the 
hoop conjecture which states that black holes with horizons form when and only when a 
mass M gets compacted into a region whose circumference C in every direction is C ^ AnM 
[5|]. This conjecture might come from the well known difference between the spherically sym- 
metric and the cylindrically symmetric systems; the spherical gravitational collapse forms 
a horizon, whereas the cylindrical gravitational collapse does not, if energy conditions on 
material fields are satisfied [sl, la]. If the hoop conjecture is correct, the following statement 
also holds: if mass M does not get compacted in some direction, there is no horizon. This 
means that if the spacetime singularity forms but its mass M does not get compacted in 
some direction, the spacetime singularity will be naked. Thus this conjecture is deeply re- 
lated also to the cosmic censorship 7!]. In order to confirm the hoop conjecture, Shapiro 
and Teukolsky performed a numerical simulation for the gravitational collapse of dust mat- 
ter with spindle-like distributions [8|, which is a relativistic counterpart of Lin-Mestel-Shu 
collapse [9[ . Their result implies that a spacetime singularity without a horizon forms from 
highly elongated distribution of matter. This is a strong candidate for a counter-example of 
the cosmic censorship conjecture. However, based on the study about an infinitesimally thin 
cylindrical shell composed of dust particles with non-vanishing angular momenta, Aposto- 
latos and Thorne argued that the effect of the rotation will halt the spindle gravitational 



collapse [10||. However, any complete answer for this issue has not been given yet [ill. Il2|. 

In this paper, we reanalyze the cylindrical shell model originally studied by Apostolatos 
and Thorne. Our study is based on the evaluation of the C-energy at the future null infinity. 



whereas the argument by Apostolatos and Thorne is on the C-energy for bounded domains. 
Due to the difference between ours and theirs, the conclusion is rather different from each 
other. 

In this paper, we adopt the geometrized units c = 1 = G, and our notation follows the 
textbook of Hawking and Ellis |13i] . 

II. CYLINDRICAL SHELL COMPOSED OF COUNTER-ROTATING PARTI- 
CLES 

A. Spacetime with whole cylinder symmetry 

In this section, we review the infinitesimally thin cylindrical shell model studied by Apos- 



tolatos and ThorneflOJ]. Hereafter, we refer to this shell as the Apostolatos-Thorne (AT)- 
shell. The AT-shell is composed of dust particles with an identical rest mass and an absolute 
value of angular momentum, and the half of those has a positive angular momentum, whereas 
the other half has negative one, such that the net angular momentum is zero, and all the 
particles remain to be on a cylindrical shell. In this case, the spacetime (M, g) has the whole 
cylinder symmetry or equivalently has the metric tensor with the local form [ij, [l5| 



g = e^^^*-'^) {-dt* ® dt* + dr* ® dr*) + r^e-^'^d^ ® d^ + e^'^dz ® dz, (1) 

where 7*, ip and r depend on t* and r*. The ranges of coordinates are given by —00 < 
t*,z < +00, < r* < +00, < (f < 27r. The axis of rotational symmetry is located at 
r* = 0. The coordinate basis d/dif and d/dz are the rotational and translational Killing 
vectors, respectively. 

If the stress-energy tensor, T, satisfies the condition 

then the Einstein equations impose the wave equation in the 2-dimensional Minkowski space- 
time on the radial coordinate function as 

(9f.r -d'^,r = 0. (3) 

The solution of the above equation takes a form 

r = fiv*)+9iw*), (4) 



where 

v* = t*+r\ w*=t*-r*, (5) 

and / and g are arbitrary functions. Here we restrict ourselves to the case that r = [const] 
hypersurface is timehke, i.e., the inequahty (9f*r)^ < {dr*r)^ holds. Then we adopt the 
metric variable, r, as a new radial coordinate, and further adopt a function defined by 

t = f{v*)-g{w*) (6) 

as the new time coordinate. The metric in this new coordinate system can be expressed as 

g = e^^"'-^^ {-dt ®dt + dr® dr) + r^e-^^d^ ® d^ + e^^dz ® dz, (7) 

where 

7 = 7*-^ln[(9.,rf-(9t.rf] (8) 

has been defined. Here only a pair of metric variables, 7 and ipy appears. 

We assume that the AT-shell is put in the vacuum spacetime. Therefore, the metric takes 
the form of Eq. (JTj) both inside and outside regions of the AT-shell. The Einstein equations 
lead to the equations for 7 and ip as 

da = 2r{dttP)dr^, (9) 

dr7 = r [{dt^y + (9,^')'] , (10) 

d^ - dl - -^^^ V- = 0. (11) 



We also assume that the space is not closed in r-direction, i.e.. 



re 



^ > for r > 0. (12) 



B. Description of the AT-shell 

The trajectory of an AT-shell in the spacetime is a timelike hypersurface. Sat- Though 
the spacetime is singular on Sat, it can consistently be treated by Israel's metric junction 
method la, [l7|, [l8|. The AT-shell, Sat, divides the spacetime into two regions. We refer to 
the inside region of Sat as M_ and the outside one as M_|_. Even if Sat is singular, we can 
require that the metric tensor g and the Killing vectors d/dip, d/dz are continuous at Sat- 



It can be easily seen that the continuity of g{d/dip, d/dip) and g{d/dz, d/dz) imphes the 
continuity of the coordinate function, r, and the metric variable, ip, across Sat- By contrast, 
the continuity of 7 across Sat is not guaranteed and this means that the coordinate function, 
t, and accordingly, the coordinate basis, d/dt, may not be continuous across Sat- Then, 
the evolution of the cylindrical AT-shell is characterized by its radial coordinate r = R{t), 
where r is the proper time naturally defined on the AT-shell. The circumferential radius of 
the AT-shell is given by TZ{t) = e~'^"*^^^i?(r), where ^s{t) is the value of ■?/' evaluated on the 
AT-shell. 

We introduce the proper reference frame of an observer riding on the AT-shell as follows, 

Eu = X±- h V—- = four velocity of the shell, (13) 

dt± or 

d d 

Eat = V— h X±^— = outward unit vector normal to the shell, (14) 

dt± or 

E^ = — ^, (15) 



r dip'' 
d_ 



where 



E- = e-*-i-, (16) 



dtA 



X± := ^ = A/e-^fe-V-.) + V^. (18) 

dr 

The subscripts -|- and — are used to denote quantities evaluated on the outer and inner faces 
of the AT-shell, respectively, if necessary. 

As mentioned above, the AT-shell is made of counter rotating dust particles which move 
along timelike geodesies whose tangents are denoted by u. By virtue of the rotational 
isometry generated by d/dip, the specific angular momentum (the angular momentum per 
unit rest mass) of each particle is conserved. Therefore, the component of u in the direction 
tangent to E^^ is given by 

g(u,E^) = ±^=:±M, (19) 

where a > is a positive constant corresponding to an absolute value of the specific angular 
momentum. Since the rest mass of each particle is conserved quantity, the shell's rest mass 
per unit Killing length is also conserved, which we denote by A. We assume that A is positive. 



i.e., A > holds. Then the surface stress-energy tensor, S, of the AT-shell is given by 

S = a(Eu(S)Eu + Y^^<P ® E J , (20) 

where we have defined 



as the surface energy density of the AT-sheU. 

In accordance with the Israel's prescription, the Einstein equations for the AT-shell reduce 
to 



K , - K_ = Stt 



S-^(TrS)/i 



(22) 



where K_|_ and K_ are the extrinsic curvatures of the AT-shell relative to the external region 
M_|_ and the internal region M_, respectively. The above equation leads to the junction 
conditions on the metric variables as 

2A 

E^^+ - E^vV-- = —-j===, (23) 

X+-X^ = -AXVl + u\ (24) 



and 



dV x.^ , ..r/x. , x2 /x. , x2i X_En^_ X_A X_X+u^ , , 



where 



E^/s = Xi^ + 1^1^ and E^/± = \/|^ + Xi^. (26) 



C. Momentarily Static and Radiation- Free Initial Data 

Here we consider the initial data with the momentarily static and radiation-free (MSRF) 
conditions 

1^ = and 5iV^ = = a>. (27) 

By solving the Einstein equations (l9l)- (iTTl) . the metric variables in M+ are given by 



7 = 7+ + fi:2ln(r/i?), (28) 

^ = 7/'i-Kln(r//2), (29) 



where 7+, ipi and k are integration constants. In contrast, the solutions in M_ should satisfy 
the regularity condition at the symmetry axis, r = 0, as 

7|r=0 = 0, dr1p\r=0 = 0. (30) 

Since ip should be continuous across the AT-shell, the metric variables in M_ are given by 

7 = 0, V' = V'i- (31) 

The junction conditions ( l23l) and ( 124|) give relationships between the integration constants 
7+, K, ipi and the quantities characterizing the AT-shell as 



7+ = -ln(l-4AVl + M2), (32) 

n = , ^"^ _ , (33) 



;i-4Av/rT^)v^ 



u 



2- 



where A := Ae '^^ is the rest mass per unit proper length. Here note that in order that 



the space is not closed in r-direction, AvT+t? < 1/4 must be satisfied. This condition 
guarantees the positivity of k, i.e., k > 0. 

As shown by Apostolatos and Thorne, Eq. ( !25l) is rewritten in the form 



where 



dV 

—— = (positive quantity) x [Acq — A] , (34) 

dr 



Ae.(«) := |YT2^- ^^^^ 

The above equation shows that if the rest mass per proper length A is greater than Aeq, 
the MSRF AT-shell starts contracting, whereas if A is smaller than Agq, it starts expanding. 
The initial data of A = Acq corresponds to the static configuration. We can easily see from 
Eq. (!33|) that k, = 2u^ in the static case. 



III. C-ENERGY ARGUMENT FOR THE FINAL CONFIGURATION 
A. C-energy 

The C-energy E has been proposed by Thorne as a quasi-local energy which is the energy 



included within the cylinder with finite radius and with unit Killing length 191]. 



-4 



^*-]^\n{{dr*rf-{dt*rf] 



(36) 



Because the external region of the AT-shell is vacuum, the expression for the C-energy 
reduces to 

E := 1 (37) 

As shown in the preceding section, the metric function, 7, diverges logarithmically in the 
limit of r — > +00 for the MSRF initial data. This means that E{t, +00) of MSRF initial 
data is necessarily infinite and thus the value of the "total" energy per unit Killing length 
is meaningless. The similar situations also appear in the Newtonian gravity; the depth of 
gravitational potential produced by cylindrically distributed matter is infinite. However, the 
temporal variation of the total energy per unit Killing length is meaningful and crucial for 
a later discussion. 

The advanced and retarded time coordinates v and w, defined by 

V = t + r, w = t — r, (38) 

are convenient to express the C-energy carried by the gravitational waves. We refer to the 
null hypersurface, v -^ +00, at infinity as the future null infinity, X"*". (About general 
arguments on the future null infinity in the spacetime with a space-translation Killing field, 
, |21j .) Due to the gravitational emissions, the C-energy will vary with w on the 



see Refs. 



20 



null hypersurface, v = [const]. The retarded time function w plays a role of a time function 
on the null hypersurface given by f = [const]. From Eqs. iQ and (TTOll . the derivative of E 
with respect to w becomes 

d^E = ^9^7 = -^{d^^y < 0. (39) 

Taking the limit of f — » +00, it can be seen that the C-energy is non- increasing function on 
X+ with respect to w, which was shown by Stachel by the argument in terms of the news 
function [3|. 

Here we study the time variation of the C-energy on X"*", which will be caused by the 
gravitational waves from the AT-shell initially in the MSRF state. Let the MSRF initial 
Cauchy surface be located at t = ti in terms both of the external and internal time coor- 
dinates, and let us refer to this Cauchy surface by Si. The initial radial coordinate of the 
AT-shell, Sat H ^1, is denoted hj R = Ri. Then the causal future of the initial location of 
the AT shell, J+(Sat H S^), is given by the condition, v^ > Vi and w_ > Wi, where Vi and 



Wi are given by 



V\ := ti + i?i, Wi := ti — Ri. 



(40) 



The portion of the spacetime, D'^{Tii) — J+(Sat H Si), remains static, while the region 
J"^(Sat n Si) will be dynamical due to the gravitational waves generated by the motion of 
the AT-shell. 

Let us describe the evolution of ip on M+ as 



^ = ^i- K\n— + 5^. 



(41) 



By virtue of the linearity of Eq. flTTl) . 6^ is subject to the same equation as that for ip as 



df - dl - -a Lu = 0. 



r I ^ip 



(42) 



Because the region D+(Si) — J+(Sat H Si) remains static, the solution in M_|_ of our interest 
takes the form 22 1 



t-r 



piO 



-A 



piO 



=^e, 



(43) 



where p{^) is a function which vanishes for C, < Wi and will decay for ^ — >■ +oo. The 
derivative of 6^ with respect to w is given by 

pio .. r pio 



0^5^ = lim- 
€^0 e 



lim - 

€^0 e 



tO+€ 



P(e + e) 



c/^ 



>v/(^-e-e)(^-0 
1 d ( piO 



d^ 



piO 



d^ 



Vi^-Oiw-0 



d^ 



dc 



Substituting the above results into Eq. (139!) . we have 



(44) 



lim ^..-E 



v-^+oo 



where 



P(w) 



1 (ip 



m^ 



(45) 



(46) 



'-00 V^^-^ d^ 
has been defined. Provided that P{w) is finite, the rate of change in E is also finite on X+. 



In order to define the quasi-local energy finite even in the limit of r ^ oo for the system 
with whole cylinder symmetry, Thorne introduced an alternative definition of the C-energy 

as 

^(new) _l(^_g-8S)_ (47) 

8 
For the MSRF initial data, this new version of the C-energy always becomes 1/8 in the limit 

of r -^ oo. For t > ti, 7 in M+ is written in the form 

7 = 7+ + fi:^ In ^ + 6^. (48) 

From Eqs. (l37|) and (l45l) . we find that 6^ at X^ is determined by 

lim d^S^ = -2P^{w), (49) 

and therefore, we have 

lim S^ = -2e{w - Wi) / P\Od^, (50) 

where 9{x) is the Heaviside's step function. From the above results, it can be seen that 6^ 



is finite. Therefore, we have 



2re 



2 



(op \ -^1^ 
1_ e-^^^++^''^d^E = 0. (51) 
V — w J 

Thus, iJ^"^^*) is constant on X"*" and this is useless for discussing how much energy is released 

into the future null infinity, X"*". For this reason, we adopt the original definition for the 

C-energy, E. 

B. Final static configuration 

The MSRF initial data with A = A; 7^ Aeq dynamically evolves. At first glance, the system 
seems to settle down into a static configuration with A = Acq by emitting or absorbing the 
gravitational waves. Thus at first, we assume that the final configuration is static and 
consider the relation between MSRF initial data and the final configuration. 

The C-energy of the initial configuration Ei and that of the final configuration E^ are 
given by 

1 



^- 4 

1 



ln(l-4Ai^/l + Mf) +'«fln-^ 



-ln( l-AA^q{ui)\/l + uf ) +fi:?ln-^ 



10 



(52) 
(53) 



Here and hereafter, quantities with a subscript "i" denote these of initial data, whereas those 
with subscript "f " denote those of the final configuration. Note that the C-energy at w = Wi 
is given hj E = E[. After the system settles down in a static configuration, the C-energy 
in the causal future of the static AT-shell is given by Eq. flS^ . Thus the difference between 
the initial C-energy at w = Wi and the final one is estimated on X+ as 

-^dE 



AE 



lim 



Wi 



dw 



dw 



lim - 

r— >+oo 4 



In 



1 - 4A,./TT^ 



[Kr 



/tf ) In r — Kj In Rf + k^ In Ri 



(54) 



l-4Aeq(Uf)v/l+W? 

We can see from the above equation that if k"^ is different from nf, the energy difference 
AE is infinite. However, we can see from Eqs. (HH!) and (!50|) that 6^ is finite and thus the 
coefficients, k, in the logarithmic term is unchanged by the gravitational emissions, i.e., 
Ki = Kf = K, hold (this was pointed out by Marder [2]). This fact is very important. On this 
ground, we can uniquely specify the static configuration realized from a given MSRF initial 
data. 

By virtue of the freedom of the constant scaling, we can assume that the initial value of 
ips, i.e., ipi, is zero, without loss of generality. Then, the initial rest mass per unit proper 
length Ai equals A. The final static configuration is characterized only by u = Uf and ips = '4'ij 
noting that the specific angular momentum a and the rest mass per unit Killing length A 
are conserved quantities. Since the final configuration is static. 



2ui 



holds. Substituting this relation into Eq. ( l33l) . we obtain 



Uf 



A 



^ (i-4Av/rT^jv/rT^ 

The final values of TZ and tps are given by using U{ in the form, 

7^f 



(55) 



(56) 



a 
ipf = In- 



A 



(57) 

Ae,(.f)- ^''^ 

Substituting Eq. (1^ into the above equations, the final value of T^f and ipf can be evaluated 

as functions of Ui. 

Here we prove the following proposition. 
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Proposition 1 Consider the MSRF initial data of the AT-shell composed of dust particles 
with non-zero specific angular momenta, i.e., a > holds. If gravitational emissions from 
the AT-shell per unit Killing length are finite, and if the initial circumferential radius TZi of 
the AT-shell is greater than the critical value given by 



7^c(wi; «, A) := W y (l - ^X^Jl + uf] ^Jl + ul (59) 

where ui = a/TZ\, i.e., TZ\ > TZc holds, then the AT-shell does not settle down into a static 
configuration. 

Proof. Suppose that the MNRS initial data of the AT-shell with the circumferential radius 
IZ = lZi settles down into a static configuration with the circumferential radius T^f . Provided 
that the gravitational emissions from the AT-shell are finite, Eqs. (1561) and (1571) imply that 
the final circumferential radius TZi is equal to IZc- Since k, remains constant, the energy 
difference ^E takes the form 



/\E = h,A + B), (60) 



where 



i-4Av/r 



u^ 



A = In ^ H , (61) 

l-4Aeq(Uf)^rT^ 

B = -2^2 In— . (62) 



Using Eqs. (l56|) . we rewrite A in the form 

4Mf A(l - A) 



A = \n 



(63) 



l + Auf{l + uf)il-\) 

where A := A/Aeq(ui). 

The condition that the space is not closed in r-direction leads the following inequality: 



1-4A-^1 + Mf >0, (64) 

which is equivalent to 

l + 4Mf(l + uf)(l-A) >0. (65) 

From Eq. (I63l) and the above inequality, A is positive if A > 1 holds. This condition, A > 1, 
is equivalent to the condition TZf <'R.i, which immediately follows from the equation 

7^? uf ^ , (l + 2wf^2 
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-(1-A), (66) 



where Eq. (156|) has been used. Hence A is positive for TZf <7l[. 
Next, from the definition of 7?. and Eq. fl58|) . we have 



From Eq. fl66|) . we have 






4(1 + nf ) + l/uf 
4(1 + nf ) + IK 
Substituting the above equation into Eq. (!67j) . we have 



^ = :)::i: 1 - (es) 



i?f 4(l + uf) + lK (1 + 2^2)2 wf^l + n; 



i?i 4{l + uf) + l/u] {l + 2ufy u^Jl+u 



f V J- -r u,f 



(69) 



It is easily seen from the above equation that Rf = Ri holds, if Uf = Ui or equivalently 
IZi = 1Z\ holds. The partial derivative of Rf/Ri with respect to U{ with Ui fixed becomes 



d ^Rf\ 4(1 + ^2) + 1/^2 i + 2Mf uf^/l + uf 



dU{\Rj [4(1+m2) + 1/m2]2 ^,6(1+ ^2)3/2 (1 + 2^2)2 

X [4(1 + u^)ut + 4(2 + 3u^)ul + 1] . (70) 

Therefore, Ri/ R\ is a decreasing function of Wf for fixed Mj, in the domain, Wi > 0. This 
implies that Ri < R\ holds, if Ui > u\ or equivalently IZi < IZi holds. Hence B is positive, if 
IZi < 1Z\ is satisfied. 

Combining the above results, it is concluded that AE > holds, i.e., the C-energy on 
X+ increases, if TZ\ > TZi = TZc holds. However this is impossible, since the C-energy is 
no n- increasing function on X+. D 

Finally, we prove the following proposition. 

Proposition 2 There exists a MSRF initial data set of an AT-shell composed of dust par- 
ticles with non-zero specific angular momenta, a > 0, which does not settle down into the 
static configuration. 

Proof. We shall show that the MSRF initial data with the condition TZ\ > TZc exists. The 
condition IZi > TZc is equivalent to 

^ > 1- (71) 



uf (l - 4Av/rT^) ^/TT^ 
We rewrite the above inequality in the form 

F{uf) > 0, (72) 

13 



where 

F{x) = X\2x + l)^-x\l+x). (73) 

In order that the space is not closed in r-direction, the inequahty 

«?<^-l (74) 

must be satisfied. (In order that this condition is satisfied, A must be smaller than 1/4.) 
Here we investigate whether Eq. (!72|) is possible in the domain < Mj < ^1/16A^ — 1. It is 
easy to see 

lim F(uf) = A^ > 0. (75) 



Since F{x) is continuous, the intersection between the domain < u^ < a/1/16A^ — 1 and 
the open neighborhood of Mi = is not empty, in which the inequality (1721) is satisfied. □ 

IV. SUMMARY AND DISCUSSION 

Apostolatos and Thorne showed that the C-energy within a bounded domain r < ro of 
MSRF configuration is minimized by the equilibrium configuration with both a and A fixed 
and further with ^(ro) = 0. They also proved that the C-energy within the domain r < tq 
of a dynamical configuration is always greater than that of a MSRF configuration that has 
the same a, A, R and ipi^ro) = 0, but different V, ifjir) and dtipir). By noticing that the 
gravitational radiation carries the C-energy from a bounded domain to outside of it, they 
inferred from the above two facts that the energy of the oscillating AT-shell is released by the 
gravitational emissions and then AT-shell will settle down into a static, equilibrium state. 
Their speculation seems to be physically reasonable, but rigorous proof has not been given 
there. 

In this paper, by assuming that gravitational radiation carries finite C-energy from the 
AT-shell to the future null infinity, we have shown that the MSRF initial data of the AT- 
shell does not necessarily settle down into the static state. Our argument is based on the 
non- increasing nature of the C-energy on the future null infinity. If the initial circumferential 
radius of the AT-shell is greater than that of the expected final static state, then the AT- 
shell can not settle down into the static state. By contrast, it is not forbidden by this 
C-energy argument that the AT-shell with an initial circumferential radius smaller than the 
expected finial static value settles down to the equilibrium configuration. We expect from 

14 



this fact that the outward motion of AT-shell is caused by gravitational emissions. If this 
expectation is real, the AT-shell with TZi > TZc will go to infinity by the secular outward 
motion. Thus the static, equilibrium configuration of the AT-shell will be unstable for the 
outward displacement. Since this secular motion may be ascribed to gravitational radiation, 
there is no Newtonian counterpart of this instability. 

At first glance, this behavior of the AT-shell is somewhat terrible, but it might not be 
so if we consider a cylinder of infinitesimally thin shell composed of dust particles but with 
finite length L which is initially much greater than its circumferential radius TZ. The motion 
of this finite cylindrical thin shell will be well approximated by the AT-shell model. Then, 
due to the instability mentioned above, TZ of the central part of the cylinder might become 
larger and larger, but finally its motion cannot be approximated by the AT-shell when TZ 
becomes comparable to L. In other words, the initial cylindrical shape might approach to the 
spherical shape by the gravitational emission. If this expectation is real, our situation seems 
to be similar to the well known phenomenon that the initial eccentric orbit of a compact star 
binary approaches to the circular orbit by the gravitational emissions [231]. However, it has 
not been rigorously shown that the above scenario is true. Rather, at present, there remains 
a possibility that the AT-shell collapses to form a naked singularity by this instability. In 
order to see whether the scenario of outward secular motion is true, we have got to investigate 
directly the solution, but this is a future work. 
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